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Abstract. In the paper we investigate a mathematical model describing the growth of tumor in the presence of im- 
mune response of a host organism. The dynamics of tumor and immune cells is based on the generic Michaelis-Menten 
kinetics depicting interaction and competition between the tumor and the immune system. The appropriate phenomeno- 
logical equation modeling cell-mediated immune surveillance against cancer is of the predator-prey form and exhibits 
bistability within a given choice of the immune response-related parameters. Under the influence of weak external 
fluctuations, the model may be analyzed in terms of a stochastic differential equation bearing the form of an over- 
damped Langevin-like dynamics in the external quasi-potential represented by a double well. We analyze properties of 
the system within the range of parameters for which the potential wells are of the same depth and when the additional 
perturbation, modeling a periodic treatment, is insufficient to overcome the barrier height and to cause cancer extinc- 
tion. In this case the presence of a small amount of noise can positively enhance the treatment, driving the system to a 
state of tumor extinction. On the other hand, however, the same noise can give rise to return effects up to a stochastic 
resonance behavior. This observation provides a quantitative analysis of mechanisms responsible for optimization of 
periodic tumor therapy in the presence of spontaneous external noise. Studying the behavior of the extinction time as a 
function of the treatment frequency, we have also found the typical resonant activation effect: For a certain frequency 
of the treatment, there exists a minimum extinction time. 

PACS. 05.40.-a-87.17.Aa-87.15.Aa 



1 Introduction 

Although cancer is a leading cause of death in the world, it 
is still little known about the mechanisms of its growth and de- 
struction. Surgery, chemo- and radio-therapies play key roles in 
treatment, but in many cases they do not represent a cure. Even 
when patients experience tumor regression, later relapse can 
occur. The need to address more successful treatment strate- 
gies is clear. Currently, efforts are made to investigate, among 
others, the methods of adoptive cellular immunotherapy Q]|2]. 
These methods of tumor treatment are based on the use of the 
injection of cultured immune cells, which have anti-tumor re- 
activity, into the tumor-bearing host. Therefore, a detailed the- 
oretical study on the mechanisms of interaction between tumor 
tissue and immune system is necessary for planning efficient 
strategies of treatment [3,4,5,6]. The immune response against 
the antigens generated by certain tumors, may be mediated by 
so called effector cells such as T-lymphocytes, macrophages 
or natural killer cells. The process of damage to tumor pro- 
ceeds via infiltration of the latter by the specialized cells which 
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subsequently develop a cytotoxic activity against the cancer 
cell-population. The series of cytotoxic reactions between the 
immune cells and the tumor tissue may be considered to be 
well approximated [3,4,5 ,6,7,8] by a saturating, enzymatic- 
like process whose time evolution equations are similar to the 
standard Michaelis-Menten kinetics: The development of tu- 
mor tissue and its reaction to immune response can be de- 
scribed in terms of a predator-prey model 16 7 8 9 10lllll[T2l 
IT^n^lT5in^rT7irT8l[Tgl[23l[24l . The population of tumor cells 
plays the role of "preys", and immune cells act as "predators". 
The activity of the predator in a certain territory, or, in this case, 
the activity of immune cells in tissue, resemble the mode of ac- 
tion of enzymes or catalysts in a chemical reaction, where the 
enzymes transform substrates in a continuous manner without 
destroying themselves. The constant immune cell population is 
assumed to act in a similar way, binding the tumor cells and 
releasing them unable to replicate. 

The network of tumor-immune system interactions is sub- 
jected to random fluctuations. The growth rate of tumor tissue 
is affected by many environmental factors, e.g. the degree of 
vascularization of tissues, the supply of oxygen, the supply of 
nutrients, the immunological state of the host, chemical agents, 
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temperature, radiations, etc. As a result of this complexity, it 
is unavoidable that in the course of time the parameters of the 
system undergo random variations which give them a stochas- 
tic character ll71[T4l[T8l[T9ll23l . 

In our previous works the effect of noise in the cancer dy- 
namics has been studied both analytically and numerically find- 
ing in its evolution the well known phenomena of resonant ac- 
tivation (RA) [25 1 and the noise enhanced stability (NES) |9] 

The aim of this study is to explore the dependence between 
the intensity of external fluctuations and the transition time 
from the state of a stable tumor to the state of its extinction 
under the influence of a moderate periodic treatment. Low in- 
tensity of therapy presents the advantage of minimizing side 
effects of treatment. The latter, if occasionally combined with 
spontaneous environmental fluctuations, can also cause the re- 
occurrence of cancer bearing a similarity to a stochastic reso- 
nance (SR) effect. Accordingly, our studies aim to understand 
a competition between the probability of the extinction and the 
reappearance of the cancer cells in the system, and to analyze 
optimal conditions of its control. 



2 The Model System 

The formulation of the model describing the growth of can- 
cerous tissue attacked by immune cytotoxic cells is based on a 
reaction scheme fl6j[7l[I3[l4j[l5][T^[l7][IS[l9l representative of 
the catalytic Michaelis-Menten scenario: 



Michaelis-Menten Potential 



X 



■2X 



X + Y 



Y + P 



(1) 



Here cancer cells X are assumed to proliferate spontaneously 
at a rate A whereas their local interactions with cytotoxic cells 

Y are modeled by a simplified kinetics with k\ standing for the 
rate of binding of immune cells to the complex Z which sub- 
sequently dissociates at a rate ku- The dissociation results in a 
product P representing dead or non-replicating tumor cells. In 
the limit case, when the production of X-type cells inhibited 
by a hyperbolic activation is the slowest process under consid- 
eration, and by assuming a conserved number of immune cells 

Y + Z = E = const, the resulting kinetics can be recast in 
the form of the first order differential equation. This can be 
interpreted as an overdamped equation of motion in a pseudo- 
potential U(x) (see M13lll4|[T5l for the details) 



dx 



dU(x) 
dx 



where 



x 2 Ox^ 

U{x) = -— + — +/3x-/3]n(x + l), 



(2) 



(3) 



x stands for the normalized molecular density of cancerous 
cells with respect to the maximum tissue capacity and the fol- 
lowing scaling relations for the kinetic parameters have been 
used 



k\x 



kiE 



t = At. 



(4) 




Fig. 1. Pseudo-potential U(x) (Eq.fJJ with parameters: j3 — 1.48, = 
0.25 that assure bistability with an almost symmetric character (equal 
depths) of both potential wells. The minima are located at Xi = and 
xs ~ 2.08, separated by a barrier at x% ~ 0.925. During the therapy 
(see the text below), this profile tilts slightly due to the driving d(t) = 
—A[l — 0(cos{2-Kvt))] which causes alternating lowering and raising 
of the barrier. The cyclic variation of the potential profile is depicted 
here for the intensity A — 0.03 and period T = 1/v = 1000 in units 
chosen accordingly to the scaling (See Eq. ©)• The minimum on the 
right side is still present during the treatment. Note a non-physical 
character of the left-hand side profile beyond xi = 0. The vertical 
draft line indicate the reflecting boundary for x — x\ in order to avoid 
this possibility on our calculations. 



Typical experimental values of the parameters are llTl fTBTTcjl : 
fci = 0.1 - 18 day- 1 , k 2 = 0.2 - 18 day" 1 , A = 0.2 - 
1.5 day' 1 . In the above deterministic model, weak fluctua- 
tions, that is temperature variations or changes in local con- 
centrations of biochemical agents, can be incorporated by in- 
cluding additional source of stochastic fluxes represented by 
an additive noise £(t) 



dx 
~dl 



dU{x) 
dx 



(5) 



which is assumed to be uncorrelated and Gaussian distributed 
with zero mean: (£(t)£(t')) = D6(t-f), and (£(*)) = 0, with 
the constraint x > 0. In the above formulation the parameter 
D 1 ' 2 stands for the intensity of the fluctuations. The resulting 
stochastic differential equation Eq.((5]l is Langevin-type. Alter- 
natively, the dynamics of the system described by the Langevin 
equation may be characterized by determining the probability 
density p(x, t) of x which carries information about the instan- 
taneous state of the system and fulfills the Fokker-Planck equa- 
tion 



d . . d 
dl pM = d-x 



dU(x) 
dx 



p(x,t) 



D d 2 P {x,t) 
2 dx 2 ' 



(6) 



The estimation of the fluctuation in cancer dynamics can be 
obtained by analyzing the cancer growth trajectories. An ex- 
ample of paper showing the time series for different patients is 
reported in (30), where we observe fluctuations in agreement 
with the range here investigated. 
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The potential profile Eq. (01 exhibits three equilibrium states 
of the system (see Fig. [TJ: 



xi = 0, 



X2 



X 3 



1-6 + y/(l + g) 2 

26 



l-6-^(TTW 
26 



(7) 
(8) 

(9) 



The essential feature captured by the model is, for a constant 
parameter 6, the /3-dependent bi-stability (GUMUlIIEIEI 
[T4lfT5l[T8lfT9l ) of stationary (long-time) solutions. In a certain 
range of values of the above parameters the model possesses 
two stable states: the state of extinction, where no tumor cells 
are present, and the state of stable tumor, where its density does 
not increase but stays at a certain constant level. The bistable 
behavior of the system can be also observed by examining the 
stationary (time-independent) solution p eg = constxe~ u ( x ^ D 
to Eq.© which, for small value of the noise intensity D 1 / 2 , is 
sharply peaked around the minima of the potential U (x). 

In the present work we are interested in situations where 
the stationarity of the system is broken due to coupling to an 
external driving representing a controllable therapy. Here, the 
influence of therapeutic factors is studied by considering a pe- 
riodic treatment (chemo- or radiation- therapy) whose action, 
at the level of the kinetic equation can be modelled by a contri- 
bution: 



d(t) = -A[l - 0(cos(2Trirt))], 



(10) 



to the RHS of Eq.[5] This term results in a periodic decrease of 
X cell concentration. The ©-symbol stands here for the Heavi- 
side function reflecting the on-off switch of the cyclic treatment 
performed with the intensity A and frequency v. Accordingly, 
the d(t) term leads to a modulation of the pseudo-potential 
whose time-dependent part is given by: 



U d {x,t) = Ax[l - 0{cos{2-Kvt))}. 



(11) 



The overall kinetic equation describing the evolution of the 
cancer cells population is then of the form 



dx 
~dt 



d{U(x) + U d {x,t)) 



dx 

= x(l - 6x) - 



(12) 



(3x 



- A[l - 0(cos(2tti4))] + 



For the purpose of these studies, we consider the situation in 
which both wells of U(x) are of equal depth, using the pa- 
rameters 6 = 0.25, and (3 = 1.48 (See Fig. IB- Since the 
cell number X cannot be negative, we investigate numerically 
the properties of the system corresponding to the Langevin 
equation ( fT2b under the constraint x > 0. In the absence of 
noise, the existence of the potential barrier prevents the trajec- 
tories (solutions to Eq. (fT2l ) from switching between potential 
wells, which correspond to well defined cancer states. If noise 
is present but there is no periodic forcing (A = 0), a sponta- 
neous tumor extinction or occurrence become possible: solu- 
tions to the Langevin equation will cross the potential barrier 
at random times whose expectation value is given by an inverse 



of Kramers rate, i.e. r oc e AU ( x ^ D where AU(x) stands for 
the height of the barrier. Under the action of periodic forcing 
—d(Ud{x, t))/dx, representing action of irradiation or chemi- 
cal therapy, the overall potential profile changes and the barrier 
will still be crossed at random times but with a preference for 
the instants of minimal barrier height. In particular, if the noise 
intensity is large enough compared to the barrier height AU, 
transitions between the wells become likely to occur twice per 
period. This phenomenon, addressed in literature [20,21] as a 
stochastic resonance (SR), is a generic phenomenon manifested 
in nonlinear systems where a weak signal can be amplified and 
optimized by the presence of noise. Obviously, if discussed in 
the context of periodic therapy administration, the SR effect 
(multiple synchronized re-crossings of the barrier) should be 
avoided as leading to possible recurrence and regrowth of tu- 
mor. 

The aim of our study is therefore to investigate the role of 
weak random contributions to deterministic dynamics of the tu- 
mor growth model in the case when the intensity (amplitude) of 
cyclical treatment is subthreshold, i.e. not sufficient to induce 
cancer extinction. A modeled action of irradiation/chemical treat- 
ment is studied based on the one-dimensional kinetic equa- 
tion (fT2l) . expressing growth of neoplastic cells population re- 
stricted to a saturation level and subject to the immunological 
surveillance mechanism of a host. Neither a particular histo- 
logical type of tumor, nor other factors determining cancer re- 
sponse to fractionated treatment (like redistribution of cells in 
various phases of cell cycle 11221 ) are specified. Within this sim- 
plified description, our main purpose is to analyze the effective- 
ness of a treatment intended to cause the extinction of a tumor 
which is already present in a tissue and whose size is described 
by the density of cancer cell population x%. 
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Fig. 2. Deterministic (D = 0) and stochastic (D > 0) trajectories of 
cancer evolution with the same parameters of Fig|T] For D = 0, the 
trajectory is localized in the cancer state at the minimum of the well. 
In contrast, a low level of noise is instead able to induce the extinction 
acts (D = 0.011). For increasing intensities of noise, the reappearance 
of cancer is detected (D = 0.016). In turn, for noise intensity close to 
the value D = 0.035 a stochastic resonance effect becomes visible. 
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3 Results 

The extinction and the reappearance of the tumor cells are stud- 
ied by the statistics of both the first extinction time (FET) and 
the first return time (FRT). We define FET as the time needed 
by a trajectory x(t) starting from the stable state x(0) — X3, to 
reach the value x — for the first time. Similarly, the FRT is 
the time needed by a trajectory x(t) to return to the initial state 
X3 for the first time, after having reached the value x = 0. To 
avoid the unphysical situation of trajectories crossing the state 
x = 0, we place there a reflecting boundary when analyzing 
the FRT, while an absorbing boundary is present in the anal- 
ysis of FET. Numerical evaluation of typical trajectories x{t) 
established for the "dose"-intensity of a treatment A = 0.03 is 
depicted in Fig. [2] Without addition of noise (D = 0) the con- 
centration of cancer cells remains close to the X3 level and the 
x{t) trajectories are localized within the right potential well (cf. 
Fig-03- in turn, a small amount of noise added to the system can 
give rise to cancer extinction whose occurrence is displayed for 
D = 0.011. 

The realistic values of the noise intensity used in this work, 
or the corresponding standard deviation a = \f~D 6 [0.005,0.32], 
can be estimated using the scaling of the Eq.|4]and the experi- 
mental values of parameters (see the previous Section), obtain- 
ing OReai ~ \Vd e [0.014,0.47] 1/day, or, in percentage of 
the maximum tumor cell density, a Real € [1.4%, 47%]. 

The statistics of such occurrences observed on the ensem- 
ble of N = 5 x 10 4 realizations x(t) is shown in Fig. [3] for 
different noise intensities. At values of D of order of 0.011 
separate peaks in histograms of FETs are visible, suggesting 
that the extinction events (passages from X3 to x = state) are 
possible only during the "treatment-on" session, beginning at 
t = T/4 up to t = 3T/4, and repeated recursively in a pe- 
riod T = 1/v = 1, 000 (in unit of 1/ A). According to the ex- 
perimental data mentioned in the previous Section, the period 
of treatment used in our simulations lies in the range (using 
T r eai = T/X [days]), T Real € [666, 5000] days. In literature 
cancer recurrence has been observed from 9 up to 17 months 
after a chemical treatment time of at least 10 months |29l . so 
the low values of that range seem to be the most realistic. 

The maximum temporal window explored is t max = 10T. 
The distinction between the two intervals ("on" and "off treat- 
ment response) becomes blurred at higher noise intensities as 
visible in Fig. [3] Moreover, by increasing the noise, the extinc- 
tion events occur at earlier times the maximum of the proba- 
bility P(t) shifts towards shorter first exit times). This means 
that the "energy supply" pumped into the system by the noise 
contributes positively to the response of the tumor cell popula- 
tion towards extinction even in the no-treatment intervals (cf. 
Fig. [5J left-bottom panel, D = 0.035). To elucidate the role 
of the spontaneous fluctuations around the deterministic tra- 
jectories, we have analyzed the statistics of cancer recurrence 
as addressed in Fig. [4] Cancer regrowth events occur first dur- 
ing the no-treatment time-interval for low level of noise in- 
tensity (D = 0.011) and, by increasing the noise-intensity, 
may be also detected during the "treatment-on" intervals. This 
apparent competition effect between the extinction and recur- 
rence events will be addressed further by analyzing the ratio 
between the mean first exit time (MFET) and the mean return 
time (MFRT). 
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Fig. 3. Probability distribution of extinction times. At very low level 
of the noise intensity, extinction events are significantly probable only 
during the treatment periods whereas no extinction is detected in the 
"treatment-off ' periods. By increasing noise intensities the extinction 
events occur also in the no-treatment intervals. Consequently, differ- 
ences between the response during "treatment-on" and "treatment-off 
intervals decrease significantly. 



At low value of the noise D = 0.011, the return probabil- 
ity is more than one order of magnitude lower than the extinc- 
tion probability estimated for a given time window (compare 
Figs. [5] and |U). This is due to the fact that the starting state for 
a temporal evolution of the X-cells concentration is X3. Con- 
sequently, the (first) returning trajectory x(t) needs more time 
to overpass the barrier, to hit the zero-concentration state and 
to continue further in a reverse direction, than an exiting trajec- 
tory which stops at x%. In addition, the potential barrier, in half 
a period, is higher in the return motion than in the exit one, in 
which the treatment facilitates the escape. On the other hand, 
at sufficiently high noise intensity, the times to cross the barrier 
in forward or backward direction are approximately equal (see 
the bottom right panel of Fig|2j. 

The total probability of the exit (P e ) and return (P r ) occur- 
rences as a function of the noise intensity calculated as the ratio 
of the number of exit/return trajectories normalized over the to- 
tal number of experiments in the observation time window, is 
shown in Fig. [5] We note that for D > 0.011 all the trajecto- 
ries leave the right well of the potential giving extinction of the 
cancer. For D > 0.024 all the trajectories experience the return 
back to the well after having reached the left minimum at x\. In 
the intermediate noise region (0.011 < D < 0.024) all the tra- 
jectories give extinction but only a certain percentage are able 
to return into the right well of the potential. The 'exit without 
return' probability P e (l — P r ) has a nonmonotonic behavior 
with an optimal value of the noise intensity (D w 0.008) not 
far from D opt shown below. The threshold lines TSe and TSr 
represent this limiting values and will be indicated as reference 
lines in the following plots. 

The above described observations are well reflected in the 
character of the mean first exit time (MFET) and the mean 
first return time (MFRT), as derived from the related proba- 
bilities P(t). Fig. [6] displays both mean times as functions of 
the noise intensity D. At very low noise intensities, an almost 
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Fig. 4. Probability density of return times. For a low level of the noise 
intensity the first return event is observed during the no-treatment time 
interval. By increasing the noise intensity the return events happens 
also during the treatment (cf.Fig[3}- At high noise intensity, the first 
return time- and the first exit time distributions tend to become similar 
each other regardless of the state of the treatment. 



deterministic character of the trajectories is mirrored in very 
high values of the MFRT and MFET, respectively. By increas- 
ing noise intensity, both MFET and MFRT decay sharply and 
MFET tends to become proportional to MFRT at a value of 
D fa 0.1. A moderate noise intensity makes the transitions be- 
tween the wells more likely to occur and thus shortens the aver- 
age extinction time of cancer. On the other hand, at higher noise 
intensities, duration of returning trajectories is (roughly) twice 
longer than duration of exiting ones. This behavior is consistent 
with the results as predicted by inspection of P(t) histograms: 
clearly, at sufficiently high noise intensities distributions of first 
exit and first return times become similar with some asymme- 
try in shape related to the asymmetry of the potential and to the 
constraints of motion imposed by the boundary condition. Con- 
sequently, the ratio MFET/MFRT is slightly higher than one 
half, as indicated in Fig|6l The two averages have been made 
using the actual number of trajectories having reached the re- 
lated boundary (x\ for the exit case, x\ and then X3 for the 
return one), excluding the trapped trajectories from the related 
calculations. The behavior of the MFET and MFRT response, 
which is different for low noise intensities, but similar for high 
noise intensities (where the two means have the same trend and 
their values tend to become proportional), accounts for the non- 
monotonic behavior of the MFET/MFRT ratio as observed in 
Fig- El In f ac t, for the chosen set of parameters, the minimum 
of MFET/MFRT is observed for D opt « 0.013 which is the 
optimum noise-intensity value giving the best compromise be- 
tween noise induced extinction of the cancer cells and their 
reappearance. 

Fig. [6] reports also the behavior of both the Kramers mean 
exit time and mean return time calculated in the case of the 
fixed potential, the expressions are the following: i) te x u = 
2Ke 2h l D / ^W'{x z )U"(x 2 )\ ; ii) r Ret = T Extt + r f with Tf = 
2ne 2h / D I y /\U"{x 1 )U"{x 2 )l where h is the height of the bar- 
rier in absence of treatment. The ratio te x u / T~R e t shows only 
a constant reference value because of the equal depth of the 
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Fig. 5. Probability of the exit and return occurrence as a function of the 
noise intensity. The 'exit without return' probability P c (l — P r ) has 
a nonmonotonic behavior with an optimal value of the noise intensity 
(D fa 0.008) not far from D opt shown in the text. 
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Fig. 6. Mean first exit time and mean first return time as a function of 
the noise intensity. The two averages have been made using the actual 
number of trajectories having reached the related boundary (xi for the 
exit case, X\ and then X3 for the return one). The ratio between both 
measures (displayed in the inset) exhibits a minimum representing the 
best stochastic compromise between effectiveness of the therapy and 
the risk of cancer recurrence. It is worth to note that in a very small 
region of noise intensity (D 6 [0.013, 0.015]) we find both the mini- 
mum value for the ratio of the MFET/MFRT and of the related width 
of the distributions (See Fig(7}- 



potential. Finally, in Figf7]the behavior of the standard devi- 
ation of the FET and FRT distributions is reported, showing 
also a well visible minimum in the ratio <je/gr at D fa 0.015, 
close to the noise intensity at which the minimum of the ra- 
tio MFET/MFRT has been observed. This means that also the 
relative precision of the exit/return times evaluations shows an 
optimal value in the same noise surroundings than the corre- 
sponding ratio of the means. Moreover, the MFET has been an- 
alyzed as a function of the frequency of the treatment, founding 
the resonant activation phenomenon ||251 (see Fig|8]l in a cer- 
tain range of the noise intensity (see also [27,28]). Note that 
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Fig. 7. Standard deviation derived from the probability densities of 
first exit and first return times. Similarly to the situation observed for 
the means, the ratio o\e I or attains a minimum within the range of in- 
vestigated noise intensities, very close to the MFET/MFRT minimum 
ztD ov t (SeeFig.[6]>. 

the numerical simulations, as discussed in this work, have been 
performed with the frequency of treatment v = 10~ 3 , corre- 
sponding to the minimum of the MFET detectable in the (D, v) 
space (Fig IS}. 



Resonant Activation Effect 

5000 - 




Fig. 8. Resonant Activation Effect due to the competitive presence 
of noise and the periodic treatment. The plot shows the mean escape 
time as a function of the noise intensity D and the frequency of the 
treatment v. The parameters are the same than those of FigQ] 



4 Conclusions 

In this paper we investigate a mathematical model describing 
the growth of tumor in the presence of the immune response 
of a host organism. The model is supplemented with periodic 



treatment and external fluctuations in the tumor growth rate. 
The formulation of the model is based on a reaction scheme |6] 
l7l [T^[T4l[T5l[T6l[T7l[T8l[T9l , representative of the catalytic Michaelis- 
Menten scenario. 

The resulting phenomenological equation modelling cell- 
mediated immune surveillance against cancer exhibits bistabil- 
ity. The two stationary points correspond to the state of a stable 
tumor and the state of its extinction. The influence of thera- 
peutic factors is introduced in the form of an on-off switch of 
the cyclic treatment, performed with a certain intensity and fre- 
quency, which (when switched on) acts in favor of a decrease 
of the tumor cell concentration. Under the influence of weak 
fluctuations, the model is analyzed in terms of a stochastic dif- 
ferential equation. It has the form of an overdamped Langevin- 
like dynamics in the external quasi-potential represented by a 
double well. We analyze properties of the system within the 
range of parameters for which the potential wells are of the 
same depth and when the periodic treatment alone (without ex- 
ternal noise) is insufficient to overcome the barrier height and 
to cause cancer extinction. 

The aim of our study was to explore the dependence between 
the intensity of fluctuations and the transition time from the sta- 
ble tumor to the extinction state. We performed a series of nu- 
merical simulations of the stochastic trajectories of the system 
and analyzed the distributions of their duration times. With- 
out addition of noise, the concentration of cancer cells remains 
close to the stable tumor state: the treatment is insufficient to 
cause the extinction. A small amount of noise added to the sys- 
tem can give rise to cancer extinction in the time intervals when 
the treatment is switched on. At higher noise intensities, the ex- 
tinction occurs even in the "treatment-off ' intervals. However, 
the presence of noise also causes the tumor recurrence. Cancer 
regrowth events occur during the no-treatment time interval for 
low noise intensity and, at higher noise intensity, may also be 
detected during the "treatment-on" intervals. At low noise in- 
tensities, the return probability can be even more than one order 
of magnitude lower than the extinction probability, but at suffi- 
ciently high noise the times to cross the barrier in both forward 
and back direction are approximately equal. 

We found that there exists an optimal noise intensity D opt 
for which the ratio MFET/MFRT of the mean first extinction 
time to the mean first return time is the smallest. Since the 
duration of particular realizations of extinction and recurrence 
processes vary in the presence of noise, we also studied the 
width of FET and FRT distributions. Interestingly enough, the 
ratio of the corresponding standard deviations (Te/ctb. a ls° nas 
a minimum at a certain noise intensity D aopt , which is close 
to D opt . Not far from the same value we also find a maximum 
for the 'extinction without return' probability as a function of 
the noise intensity. Our study shows that the presence of noise 
gives the possibility of tumor extinction even at a weak radi- 
ation or chemical treatment. Treatment strategies which take 
into account the noise contribution can reduce the radiation 
(or chemotherapeutic) dose. On the other hand, the presence 
of noise implies also the possibility of a recurrence. We show, 
however, that it is possible to observe an optimal noise value in 
the competition between the extinction and return events. 
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